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{J*} ' Abstract 

We have reviewed recent developments of some aspects of optical spatial solitons in photorefrac- 
tive media. Underlying principles governing the dynamics of photorefractive nonlinearity have been 
discussed using band transport model. Nonlinear dynamical equations for propagating solitons have 
O ■ been derived considering single as well as two-photon photorefractive processes. Fundamental proper- 

ties of three types of solitons, particularly, screening, photovoltaic and screening photovoltaic solitons 
^ ' have been considered. For each type of solitons, three different configurations i.e., bright, dark and 

■ gray varieties have been considered. Mechanisms of formation of these solitons due to single as well as 

O , two-photon photorefractive processes have been considered and their self bending discussed. Vector 

• i-H 

solitons, particularly, incoherently coupled solitons due to single photon and two-photon photorefrac- 
tive phenomena have been highlighted. Existence of some missing solitons have been also pointed 
O t out. 
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1 Introduction 



The advent of nonlinear optics has paved the way to a number of fundamental discoveries. Exotic 
one, among many of these, is the discovery of optical solitons flUl^. These are optical pulses or beams 
which are able to propagate without broadening and distortion. Optical solitons are envelope solitons, 
and the term soliton was first coined by Zabusky and Kruskal in 1965 to reflect the particle like nature 
of these waves that remain intact even after mutual collision Q . Optical solitons are extensively studied 
topic not solely due to their mathematical and physical elegance but as well as due to the possibility 
of real life applications. They have been contemplated as the building blocks of soliton based optical 
communication systems, signal processing, all optical switching, all optical devices etc. 



In optics, three different types of solitons are known till date, these are, temporal [2|, |3j, spa- 
tial 0, @] and spatio-temporal solitons 0, @] . Temporal solitons are short optical pulses which maintain 
their temporal shape while propagating over long distance. The way a temporal optical soliton is estab- 
lished is that a nonlinear pulse sets out in dispersive medium and develops a chirp. Then the dispersion 
produces a chirp of opposite sign. A temporal soliton pulse results due to the balancing of these opposite 
chirps across the width of the pulse, which arise from the material dispersion and nonlinearity. These 
opposite chirps balance each other when dispersion is completely canceled by the nonlinearity of the 
medium. Temporal solitons are routinely generated in optical fibers [2| and they are backbone of soliton 
based optical communication systems, soliton lasers etc. In contrast, optical spatial solitons are beams of 
electromagnetic energy that rely upon balancing diffraction and nonlinearity to retain their shape. While 
propagating in the nonlinear medium, the optical beam modifies the refractive index of the medium in 
such a way that the spreading due to diffraction is eliminated. Thus, optical beam induces a nonlinear 
waveguide and at the same time guided in the waveguide it has induced. This means soliton is a guided 
mode of the nonlinear waveguide induced by it. 

Though temporal solitons can be easily generated in optical fibers, generation of spatial optical soli- 
tons is a much more difficult task. For example, in silica glass, the nonlinearity is proportional to light 
intensity and the value of the proportionality constant n-i is of the order of 10 -16 cm 2 /W only. Therefore, 
in order to compensate for the beam spreading due to diffraction, which is a large effect, required optical 
nonlinearity is very large, and, consequently optical power density is also large Another impediment in 
detecting spatial solitons, in bulk Kerr nonlinear media, is the catastrophic collapse of the optical beam, 
which is inevitable with Kerr nonlinearity. Discovery of non Kerr nonlinearities, whose mechanism is 
different from Kerr nonlinearity, has lead to the revelation of stable three dimensional soliton formations 
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without catastrophic collapse. These nonlinearities are photorefractive nonlinearity [Tol. Il8l ] , quadratic 

With the iden- 



nonlinearity 19h22|| and resonant electronic nonlinearity in atoms or molecules 
tification of photorefractive nonlinearity, which possesses strong nonlinear optical response, it is possible 
to create optical solitons at very low light intensity. Spatial photorefractive optical solitons possess some 
unique properties which make them attractive in several applications, such as, all optical switching and 



routing, interconnects, parallel computing, optical storage etc [15l Il6l 1281. |35|. They are also promising 
for experimental verification of theoretical models, since, they can be created at very low power. In the 
present article, we confine our discussion on the properties of photorefractive spatial solitons. 
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2 Photorefractive Effect 



The photorefractive (PR) effect is the change in refractive index of certain electro-optic materials owing 
to optically induced redistribution of charge carriers. Light induced refractive index change occurs owing 
to the creation of space charge field, which is formed due to nonuniform light intensity. Originally the 
photorefractive effect was considered to be undesirable, since this leads to scattering and distortion of 
collimated optical beams [36J. Soon it was realized that these materials have potential applications in 



holography [33|, optical phase conjugation [32|, optical signal processing and optical storage [34j, 1351 ] . 
Photorefractive materials are classified in three different categories. Most commonly used photore- 
fractive materials are inorganics, such as, LiNbO^, BaTiOz, Sr x Bri^ x Nb 2 OQ, KNbO^, Bii 2 Si02o, 
Bi 12 Ti02o(BTO) etc. Semiconductors have large carrier mobility that produces fast dielectric response 
time, which is important for fast image processing. Therefore, photorefractive semiconductors, such as, 
GaAs, InP, CdTe etc., complement the photorefractive ferroelectrics with the potential of fast holographic 



processing of optical information. Polymers also show strong photorefractive effect [37H40l |. They are 
easy to produce and PR effect in polymers appear only if a high voltage is applied. Strong photorefrac- 
tive pattern can be erased easily in polymers by decreasing the applied voltage. Polymers show good 
temperature stability, and for a given applied voltage, they usually show stronger refractive index change 
in comparison to inorganic crystals with equal doping densities. 



3 Origin of Photorefractive Nonlinearity 



In a photorefractive material, the spatial distribution of intensity of the optical field gives rise to an 
inhomogeneous excitation of charge carriers. These charge carriers migrate due to drift and or diffusion 
and produce a space ch arg e field, whose associated electric field modifies the refractive index of the crystal 



32l [3J, |35| . For a noncentrosymmetric photorefractive crystal, the refractive index 



via Pockel's effect 
change due to the linear electro-optic effect is given by 32|. |34| 



1 „ 

An = --n r eff E sc , 



(1) 



where n is the average refractive index, r e f / is the effective linear electro-optic coefficient which depends 
on the orientation of the crystal and polarization of light, and E sc is the space charge field. A unique 
feature of photorefractive materials is their ability to exhibit both self focusing and defocusing nonlinearity 
in the same crystal. This is achieved by changing the polarity of the biased field, which in turn changes 
the polarity of the space charge field E sc . Hence, the same crystal can be used to generate either 
bright ( require self focusing nonlinearity) or dark and gray solitons ( require defocusing nonlinearity). 
Photorefractive nonlinearity is also wavelength sensitive, thus, it is possible to generate solitons at one 
wavelength and then use the soliton supported channel to guide another beam at different wavelength. 



4 Band Transport Model 

The most widely accepted theoretical formulation of photorefractive phenomenon is described by Kukhtarev- 
Vinetskii band transport model [4l|. A schematic diagram of this model is shown in figure (1), where a 
PR crystal is being illuminated by an optical beam with nonuniform intensity. 

Insert Figure (1) here 

The photorefractive medium, at the ground level, has completely full valance band and an empty con- 
duction band. The material has both donor and acceptor centers, uniformly distributed, whose energy 
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states lies somewhere in the middle of the band gap. The acceptor electron states are at much lower 
energy in comparison to that of the donor electron states. The presence of a nonuniform light beam 
excites unionized donor impurities, creates charge carriers which move to the conduction band where 
they are free to move to diffuse or to drift under the combined influence of self generated and external 
electric field and are finally trapped by the acceptors. During this process, some of the electrons are 
captured by ionized donors and thus are neutralized. In the steady state, the process leads to the charge 
separation, which tends to be positive in the illuminated region and vanish in the dark region. 

We assume that the donor impurity density is No and acceptor impurity density Na- If the 
ionized donor density is N^, then the rate of electron generation due to light and thermal processess is 
(sjJ + Pt)(Nd — Np), where Sj is the cross section of photoexcitation, / is the intensity of light which 
is written in terms of Poynting flux / = ^l*!"! 2 , $ is the electric field of light, eo is the free space 
permittivity, and (3t is the rate of thermal generation. If TV is the electron density and 7^ the electron 
trap recombination coefficient, then, the rate of recombination of ionized donars with free electrons is 
ryRNNp. Thus, the rate equation describing the donar ionization is given by 

8N + 

-^L = ( Sl I + l 3 T )(N D -N+)- lR NN+. (2) 

The electron concentration is affected by recombination with ionized donors and due to migration 
of electrons, resulting in an electron current with current density J, hence, electron continuity equation 
turns out to be 

and 

~1f = eNn'i + k B T(j^N + k pSi (N D - N+) Ic, (4) 

where the current density J is the sum of the contributions from the drift, diffusion and photovoltaic 
effect; e is the electronic charge, /x is the electron mobility, ks is the Boltzmann constant, T is electron 
temperature, k p is the photovoltaic constant and c is the unit vector in the direction of c-axis of the PR 
crystal. E is the total electric field including the one externally applied and that associated with the 
generated space charge. The redistribution of electrical charges and the creation of space charge field 
obey the Poisson's equation, therefore, 

^ • = p, (5) 

and the charge density p is given by 

p{l?) = e(N+-N A -N). (6) 

Equations (2) -(6) can be solved to find out the space charge field E sc and subsequently the optical 
nonlinearity in the photorefractive media. 



5 Space Charge Field 

To estimate the nonlinear index change in photorefractive media due to the presence of nonuniform optical 
field, we need to calculate the screening electric field E sc . The response of a photorefractive material to 
the applied optical field is anisotropic and it is nonlocal function of light intensity. Anisotropy does not 
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allow radially symmetric photorefractive solitons 42|, |43[ . To formulate a simple problem, and since, most 



of the experimental investigations on photorefractive solitons are one dimensional waves, it is appropriate 
to find material response in one dimension ( say x only). Steady state photorefractive solitons may be 
obtained under homogeneous background illumination, which enhances dark conductivity of the crystal. 
In the steady state, induced space charge field E sc can be obtained from the set of rate, continuity, current 



equations and Gauss law. In the steady state, and in one dimension, these equations reduce to 44-3 



Si (I + I d )(N D -N+)-^ R NN+=0, (7) 

J = eNfiE sc + k B Tfi^- + k pSl (N D - N+) I, (9) 

where e r is the relative static permittivity; Id{= — ) is the so called dark irradiance that phenomeno- 
logically accounts for the rate of thermally generated electrons. This is also the homogeneous intensity 
that controls the conductivity of the crystal. In most cases, the optical intensity / is such that for electron 
dominated photo-refraction, N <C Np, N -C Na and Na <C N^. Under this usually valid situation, the 
space charge field E sc is related to the optical intensity / through 

F - F i °° + Id i f J °°~ 7 kBT 1 dI nil 

where E a is the external bias field to the photorefractive crystal, E p — kpl ^ A is the photovoltaic 
field. In addition, we have assumed that the power density of the optical field attains asymptotically a 
constant value at x — > ±oo i.e., I(x — > ±oo, z) — loo- Moreover, in the region of constant illumination, 
equations(7)-(10) require that the space charge field E sc is independent of x i.e., E sc (x —¥ ±oo, z) = E a . 



6 Photorefractive Nonlinearity 

The space charge induced change An in the refractive index n is obtained as 



|An| = ~n 3 r ef f 



p loo + Id , p loo -I k B T 1 dl 
I + I D + p I + Id e I + I D dx 



(12) 



It is evident from above expression that the change in refractive index is intensity dependent, i.e., under 
the action of nonuniform illumination the photorefractive crystal has become an optical nonlinear medium. 
The first term in the above equation is generally known as the screening term. In nonphotovoltaic 
photorefractive material, this is the main term which is responsible for soliton formation. The space 
charge redistribution in a photorefractive crystal is caused mainly by the drift of photoexcited charges 
under a biasing electric field. This mechanism leads directly to a local change of refractive index and is 
responsible for self-focusing of optical beams. The second term is due to photovoltaic effect, which leads to 
the formation of photovoltaic solitons. First two terms do not involve spatial integration, therefore, they 
are local nonlinear terms. Functional form of both terms are different from Kerr nonlinearity. These terms 
are similar to that of saturable nonlinearity and explain why the collapse phenomenon is not observable 
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in photorefractive materials. Besides drift mechanism, transport of charge carriers occurs also due to 
diffusion process. This process results in the nonlocal contribution to the refractive index change. The 
last term is due to diffusion. The strength of the diffusion effect is determined by the width of the soliton 
forming optical beam. In case of strong bias field and relatively wide beams, the diffusion term is often 
neglected. However, its contribution can become significant for very narrow spatial solitons or optical 
beams. When diffusion is significant, it is responsible for deflection of photorefractive solitons 44 



7 Classification of Photorefractive Solitons 

7.1 Solitons Due to Single Photon Photorefractive Phenomenon 

Till date, three different types of steady state spatial solitons have been predicted in photorefractive 
media, which owe their existence due to single photon photorefractive phenomenon. Photorefractive 
screening solitons were identified first. In the steady state, both bright and dark screening solitons (SS) 
are possible when an external bias voltage is appropriately applied to a non-photovoltaic photorefractive 



crystal 47-5(3]. Intuitively, the formation of bright screening photorefractive solitons can be understood 



as follows. When a narrow optical beam propagates through a biased photorefractive crystal, as a result 
of illumination, conductivity in the illuminated region increases and the resistivity decreases. Since the 
resistivity is not uniform across the crystal, the voltage drops primarily across the non illuminated region 
and voltage drop is minimum in the high intensity region. This leads to the formation of large space charge 
field in the dark region and much lower field in the illuminated region. The applied field is thus partially 
screened by the space charges induced by the soliton forming optical beam. The refractive index changes 
which is proportional to this space charge field. The balance of self diffraction of the optical beam by the 
focusing effect of the space charge field induced nonlinearity leads to the formation of spatial solitons. 
These screening solitons were first predicted by Segev et. al. ^3]> whereas, the experimental observation 
of bright SS were reported by M.Shih et. al. [48j and those of dark SS were reported by Z. Chen et. al. 



The second type PR soliton is the photovoltaic solitons [57H6l|, the formation of which, however, requires 
an unbiased PR crystal that exhibits photovoltaic effect, i.e., generation of dc current in a medium illu- 
minated by a light beam. The photovoltaic (PV) solitons owe their existence to bulk photovoltaic effect, 
which creates the space charge field, that, in turn modifies the refractive index and gives rise to spatial 
solitons. These solitons were first predicted by G. C. Valey et. al. 57] and observed experimentally in 



ID by M.Taya et.al. [58| and in 2D by Z.Chen et. al. [59[. Two dimensional bright photovoltaic spatial 
solitons were also observed in a Cu:KNSBN crystal by She et. al. [(ij . The observed spatial solitons 
were broader than those predicted by considering E sc due to signal beam alone. This was satisfactorily 
explained by an equivalent field induced by the background field. 

The third type of photorefractive solitons arises when an electric field is applied to a photovoltaic pho- 



torefractive crystal [3l|, |60|, |6l|, |63|, |64| . These solitons, owe their existence to both photovoltaic effect and 
spatially nonuniform screening of the applied field, and, are also known as screening photovoltaic (SP) 
solitons. It has been verified that, if the bias field is much stronger than the photovoltaic field, then, the 
screening photovoltaic solitons are just like screening solitons. On the other hand, if the applied field is 
absent, then SP solitons degenerate into photovoltaic solitons in the closed circuit condition. The first 
observation of bright photovoltaic screening solitons in LiNbOy, was reported by E. Fazio et. al. [3lj j . 
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7.2 Nonlinear Equation for Solitons 

In order to develop a semi analytical theory, one dimensional reduction is introduced in the subsequent 
discussion. The optical beam is such that no y dynamics is involved and it is permitted to diffract only 
along x direction. Electric field E of the optical wave and the bias field Eq are directed along the x 
axis which is also the direction of the crystalline c-axis. Under this special arrangement, the perturbed 
extraordinary refractive index h e is given by 
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46, 



(h e ) = n e - n e r eff E sc , (13) 

where n e is the unperturbed extraordinary index of refraction. This arrangement allows us to describe 
the optical beam propagation using Helmholtz's equation for the electric field i^, which is as follows: 

V 2 E* + {k h e ) 2 ~^ = 0, (14) 

where ko = jL anc l \ Q j s the free space wavelength of the optical field. Moreover, if we assume 
z)exp[i(kz — uit)], where k = kon e and employ slowly varying envelope approximation for 
the amplitude $, then, the Helmholtz's equation can be reduced to the following parabolic equation: 

<9$ 1 <9 2 $ 1 
i- 



■8z ' 2k 6x* 2^/^ = 0. (15) 
The space charge field that has been evaluated earlier in section(5) can be employed in the above equation 
to obtain following nonlinear Schrodinger equation: 

dA ld 2 A . A (p-\A\ 2 )A r A d\A\ 2 , . 

l a[ + 2a?- p{1 + p) T+W a i + \a\ 2 + s i + \a\ 2 ds = °' (16) 

where C = fe^ij> s = ^. P = 77' P = (kox Q ) 2 {n 4 e r eff /2)E , a = {k x ) 2 (n 4 e r eff /2)E p , A = 
$ and S — (k 2 xor e ffn^kBT)/(2e). Equation (16) gives rise to varieties of solitons depending on 



specific experimental situations. Important parameters which classify these solitons and govern their 
dynamics are a and /?. The parameter 5, which is associated with the diffusion term, is not directly re- 
sponsible for soliton formation. The diffusion processes is primarily responsible for bending of trajectories 
of propagating solitons, hence, large value of S influences the trajectory of bending. The dimensionless 
parameter j3 can be positive or negative depending on the sign of Eq i.e., the polarity of external applied 
field. In non-photovoltaic photorefractive media a = 0, therefore, if we neglect the diffusion term, which 
is small indeed, and introduce only first order correction, then /3 is the parameter which governs the 
soliton formation. The relevant dynamical equation in non-photovoltaic photorefractive media turns out 
to be 

.dA Id 2 A OM t . A 

^ + 2^-^ + p) TTW = °- (17) 

This is the modified nonlinear Schrodinger equation(MNLSE) with saturable nonlinearity which is 
not integrable. The term 1/(1 + \A\ 2 ) represents local saturable nonlinear change of refractive index of 
the crystal induced by the optical beam. The saturating nature of the nonlinearity will be more clearly 
evident if we make the transformation A — u exp[—i/3(l + p)£], in which case above equation reduces to 

du 1 d 2 u „, , lul 2 , , 

^ + 2d^ + ^ + ^TTW u = - (18) 
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It will be shown subsequently that a positive /3 is essential for creation of bright spatial solitons, 
whereas, a negative /3 leads to the formation of dark or gray solitons. Thus, by changing the polarity of 
the external applied field, it is possible to create bright as well as dark solitons in the same photorefractive 
media. Equation(17) has been extensively investigated for bright, dark as well as gray solitons 



60 



4G 



6l| . In the next section, we consider bright screening spatial solitons using a method outlin ed by 



Christodoulides and Carvalho 



111 



137 



145 



146|. 



and subsequently employed in a large number of investigations 
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7.3 Bright Screening Spatial Solitons 

We fist consider bright screening solitons for which the soliton forming beam should vanish at infinity 
(s — > ±oo), and thus, = p = 0. Therefore, bright type solitary waves should satisfy 

dA 1 d 2 A A 

^ + 2^-^TTW = °- (19) 

We can obtain stationary bright solitary wave solutions by expressing the beam envelope as A(s, £) = 
piy(s) exp(zi/£), where v is the nonlinear shift in propagation constant and y(s) is a normalized real 
function which is bounded as, < y(s) < 1. The quantity p represents the ratio of the peak intensity 
(I m ax) to the dark irradiance Id, where I m ax = I(s = 0). Furthermore, for bright solitons, we require 
y(0) = 1, y(0) = and y(s — > ±oo) = . Substitution of the ansatz for A(s,£) into equation(19) yields 

-l-2vy- 2/3—^ = 0. (20) 
ds z 1 + py z 

Integration of above equation once and use of boundary condition yields 

i/ = -(/9/p)lii(l+p), (21) 

and 

=(2/3/p)[ln(l+ ra 2 )-y 2 ln(l+p)]. (22) 

By virtue of integration of above equation, we immediately obtain 

1 f 1 p 1/2 
S = ± j2^J y [ln(l+pP)-y 2 ln(l+p)]i/2 dy - (23) 

The nature of above integrand is such that it does not provide any closed form solution. Nevertheless, 
the normalized bright profile y(s) can be easily obtained by the use of simple numerical procedure. It 
can be easily shown that the quantity in the square bracket in equation(23) is always positive for all 
values of y 2 between < y(s) < 1. Therefore, the bright screening photorefractive solitons can exist in 
a medium only when /3 > i.e., Eq is positive. For a given value of /3, the functional form of y can be 
obtained for different p which determines soliton profile. For a given physical system, the spatial beam 
width of these solitons depends on two parameters Eq and p. For illustration, we take SBN crystal with 
following parameters n e = 2.35, = 224 x 10 -12 m/V. Operating wavelength Ao = 0.5fim, xq — 20fim 
and £o = 2x 10 5 V/m. With these parameters, value of f3 ~ 43. Figure (2) depicts typical normalized 
intensity profiles of bright solitons. 
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Insert Figure (2) here 



7.3.1 Bistable Screening Solitons 

Equation(19) possesses several conserved quantities, one such quantity is P = f \A\ 2 ds, which can be 
identified as the total power of the soliton forming optical beam. Numerically evaluated soliton profile 
\A\ 2 = P\y( s )\ 2 can be employed to calculate P. Solitons obtained from equation(23) are stable since 
they obey Vakhitiov and Kolokolov j^] stability criteria i.e., dPjdv > 0. These soliton profiles can be 
also employed to find out spatial width tfwhm ( full width at half maximum) of solitons. Figure (3) 
demonstrates the variation of spatial width tfwhm with p. This figure signifies the existence of two- 
state solitons, also known as bistable solitons, i.e., two solitons possessin g sa me spatial width but different 



power. Similar bistable solitons were earlier predicted in doped fibers 66, 67]. However, these bistable 



solitons are different from those which were predicted by Kaplan and others [68[ , where two solitons with 
same power possessing two different nonlinear propagation constant. 

Insert Figure (3) here 

Another point worth mentioning is that, the tfwhm vs p curve in figure(3) possesses local minimum, 
hence, it is evident that these solitons can exist only if their spatial width is above certain minimum 
value. This minimum value increases with the decrease in the value of (3. For illustration, the shapes 
of a typical pair of bistable solitons with same tfwh m but different peak power have been depicted in 
figure(4). The dynamical behavior of these bistable solitons, while propagating, can be examined by full 
numerical simulation of equation(19), which confirms their stability. 



Insert Figure (4) here 



7.4 Dark Screening Solitons 



Equation (17) also yields dark solitary wave solutions [46|,|69(, which exhibit anisotropic field profiles with 
respect to s. These solitons are embedded in a constant intensity background, therefore, 1^ is finite, 
hence, p is also finite. To obtain stationary solutions, we assume A(s,£) = p 1 ' 2 y(s)exp(iv^), where, like 
earlier case, v is the nonlinear shift in propagation constant and y(s) is a normalized real odd function 
of s. The profile y(s) should satisfy following properties: y(0) — 0,y(s — > ±oo) = ±1, = — as 
s — > ±oo. Substituting the form of A in equation(17) we obtain, 

±]L - 2vy- 2(3(1 + p)— ^ = 0. (24) 
ds z 1 + py^ 

Integrating above equation once and employing boundary condition, we get 

v = -fi. (25) 



Following similar procedure as employed earlier, we immediately obtain 

= | i f° € 

(-2/3)V2 J y w _ l) _(i±P) ln i±ml ]1/ 2- 1 > 

The quantity within the square bracket in equation(26) is always positive for all values of y 2 < 1, thus 
i.e., Eq must be negative so that r.h.s of equation(26) is not imaginary. An important point to note 
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is that, for a particular type of photorefractive material, for example SBN, if positive polarity of Eq is 
required for bright solitons, then one can observe dark solitons by changing polarity of Eq . Spatial width 
of these solitons depends on only two variables j3 and p i.e., Eq and 1^. It has been confirmed that, 
unlike their bright counterpart, dark screening photovoltaic solitons do not possess bistable property. For 
illustration, we take the same SBN crystal with other parameters unchanged, except in the present case 
Eq = —2 x 10 5 V/m. Therefore f3 sa —43. Figure (5) depicts normalized intensity profiles of dark solitons 
which are numerically identified using equation(26). 

Insert Figure (5) here 



7.5 Gray Screening Solitons 



Besides bright and dark solitons, equation (17) also admits another interesting class of solitary waves, 
which are known as gray solitons 46j. In this case too, wave power density attains a constant value at 
infinity i.e., 1^ is finite, and hence, p is finite. To obtain stationary solutions, we assume 



A(s,$) = p 1/2 y(s) exp[«K 



Jds 
y 2 (s) 



(27) 



where v is again the nonlinear shift in propagation constant, y(s) is a normalized real even function 
of s and J is a real constant to be determined. The normalized real function satisfies the boundary 
conditions y 2 (0) = m(0 < to < 1), i.e., the intensity is finite at the origin, y(0) = 0, y(s — > ±oo) = 1 and 
all derivatives of y(s) are zero at infinity. The parameter m describes grayness, i.e., the intensity 1(0) at 
the beam center is 1(0) = m/oo. Substitution of the above ansatz for A in equation(17) yeilds 



-2vy -2/3(1 + p) T 



py 



J 2 
y 3 



o. 



(28) 



Employing boundary conditions on y at infinity we obtain 

J 2 = —2(v + f3). (29) 
Integrating equation(28) once and employing appropriate boundary condition, we immediately obtain 

(30) 



(to - If 



m( l + P) , fl + pm 

In I — | + 1 — to 



1 + p 



Finally 



(31) 



The normalized amplitude y(s) can be obtained by numerical integration of above equation. Note 
that dark solitons are a generalization of these gray solitons. Unlike bright or dark solitons, the phase 
of gray solitons is not constant across s, instead varies across s. Existence of these solitary waves are 
possible only when /3 < 1 and to < 1. 
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7.6 Self Deflection of Bright Screening Solitons 

In the foregoing discussion we have neglected diffusion, however, the effect of diffusion cannot be ne- 
glected when solitons spatial width is comparable with the diffusion length. The diffusion process in- 
troduces asymmetric contribution in the refractive index change which causes solitons to deflect during 
propagation. Results of large number of investigations addressing the deflection of photorefractive spatial 
solitons in both non-centrosymmetric and centrosymmetric photorefractive crystals are now available 
l^ ]. Several authors have investigated self bending phenomenon using perturbative procedure 
In this section, we employ a method of nonlinear optics 85, which is different from perturbative 




approach. To begin with, we take finite 6 and use the following equation to study the self bending of 
screening bright spatial solitons 



.d A Id 2 A 



A 



+ 8- 



A d(\A\ 2 ) 



d£ 2 ds 2 ^ {l + \A\ 2 ) {l + \A\ 2 ) 
To obtain stationary solitary waves, we make use of the ansatz 



ds 



= o. 



(32) 



A( s ,£)=A (s,£)expHQ( s ,£)], 
in equation(32). A straightforward calculation yields following equations: 

dA Q dA Q on 1 d 2 n 

:Ao^r^r — 0, 



<9£ ds ds 2 u ds 2 



(33) 



(34) 



and 



dQ ld 2 A 1 fdn^ 2 



A 



S d(A 2 ) 

1 + A 2 ^ 1 + A 2 ds 



(35) 



We look for a self-similar solution of (34) and (35) of the form 



A 



Mm 



exp 



21 



2r 2 f 2 (0 



nM = (« s °(0) 2 Ai(g) + {s So(€ ))A 2 (0 + A 3 (£), 



(36) 
(37) 



Ai = 



1 df 



and 



A - ds ° 
d£ 



(38) 



where, ro is a constant and /(£) is a parameter which together with tq describe spatial width; in 
particular, ro/(£) is the spatial width of the soliton and A3 is an arbitrary longitudinal phase function. 
So(£) is the location of the center of the soliton. For a nondiverging/ nonconverging soliton, /(£) = 1. 
Moreover, we assume that initially solitons are nondiverging i.e., ^jp- = at £ = 0. Substituting for Aq 
and fl in equation(35) and equating coefficients of s and s 2 on both sides, we obtain 



de 



rW r 2 f 



9 A 2 

„2 f 2 



and 



1 



A 2 



(39) 



d 2 s (0 = _ § 2A 2 



r 2 of 3 



A 2 



(40) 
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Equation(39) describes the dynamics of the width of soliton as it propagates in the medium, while 
equation(40) governs the dynamics of the centre of the soliton. In order to find out how a stationary 
soliton deviates from its initial propagation direction, we first solve equation(39) for stationary soliton 
states. From equation(39), condition for stationary soliton states can be obtained as 



(1 



2/L4 2 



no 



1/2 



(41) 



Insert Figure ( 6) here 

Figure (6) depicts variation of solitons spatial width with the normalized peak power . The curve 
in figure (6) is the existence curve of stationary solitons. Each point on this curve signifies the existence of 
a stationary soliton with a given spatial width and peak power. A stationary soliton of specific power and 
width as described by equation(41) deviates from its initial path which can be found out by integrating 
equation(40). The equation of trajectory of the center of soliton is 

so(£) = -|£ 2 + soo, (42) 

where 8 = S^r-Cl+Ann), sqq — so(£ = 0), moreover / = 1, since we are only interested in stationary 
solitons whose spatial shape remain invariant. Thus, the beam center follows a parabolic trajectory. 
The displacement of the soliton center with propagation distance has been depicted in figure (7). It is 
evident that the peak power of soliton influences lateral displacement of the soliton centre. The lateral 
displacement suffered by the spatial soliton is given by ^£ 2 - Equation(42) implies that the angular 
displacement of the soliton center shifts linearly with the propagation distance £. The more explicit 
expression of the angular displacement, i.e., the angle between the center of the soliton and the z axis 
can be easily obtained as 0£. 

Insert Figure (7) here 

7.7 Photovoltaic and Screening Photovoltaic Solitons 

Steady state photovoltaic solitons can be created in a photovoltaic photorefractive crystal without a bias 
field. These photovoltaic solitons result from the photovoltaic effect 57 - 6l| • Recently, it has been 



predicted theoretically that the screening-photovoltaic(SP) solitons are observable in the steady state 



when an external electric field is applied to a photovoltaic photorefractive crystal [3l|, [78J. These SP 
solitons result from both the photovoltaic effect and spatially nonuniform screening of the applied field. 

If the bias field is much stronger in comparison to the photovoltaic field, then the SP solitons are just 
like the screening solitons. In absence of the applied field, the SP solitons degenerate into photovoltaic 
solitons in the close-circuit condition. In other words, a closed-circuit photovoltaic soliton or screening 
soliton is a special case of the SP soliton. Thus, in the subsequent analysis, we develop theory for 
SP solitons and obtain solutions of photovoltaic solitons as a special case. When the diffusion process 
is ignored, the dynamics of these steady-state screening PV solitons ( bright, dark and gray) can be 
examined 



60 



78l l88l | using the following equation 



.OA Id 2 A „.„ , A (p-\A\ 2 )A n 

' ac + 2 a* ~ p{1 + "> (T+m - a jrrwy = °- (43) 
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Adopting; the procedure, employed earlier, the bright soliton profile of screening photovoltaic soli- 
tons [60( turns out to be 

s(2(q + P)) 1 ' 2 =± ( dy, (44) 

J y [ln(l + py 2 ) - y 2 ln(l + p)\ L ' 2 

The quantity within the bracket in the right hand side is always positive for y(s) < 1, thus, (a + /?) 
must be positive for bright screening photovoltaic solitons. From equation(44), it is evident that, these 
bright SP solitons result both from the photovoltaic effect (a ^ 0) and from spatially nonuniform screen- 
ing (/3 ^ 0) of the applied electric field in a biased photovoltaic-photorefractive crystal. Formation of 
these solitons depends not only on the external bias field but also on the photovoltaic field. When we 
set a = 0, these solitons are just like screening solitons in a biased nonphotovoltaic photorefractive crys- 
tal 40| . In addition, when /3 = , we obtain expression of bright photovoltaic solitons in the close circuit 



realization 



611 ] . Thus, these SP solitons differ both from screening solitons in a biased nonphotovoltaic 
photorefractive crystal and from photovoltaic solitons in a photovoltaic photorefractive crystal without an 
external bias field. One important point to note is that the experimental conditions are different for cre- 
ation of screening, PV and SP solitons. SP solitons can be created in a biased photovoltaic-photorefractive 
crystal, whereas, creation of screening solitons are possible in biased nonphotovoltaic-photorefractive crys- 
tals. The PV solitons can be created in photovoltaic-photorefractive crystal without an external bias field. 

We can also derive dark solitons from equation (43), the normalized dark field profile can be easily 
obtained as 

[-2(a + P)]W 8 = ±[° ^ -TTj, (45) 



6P-l)-±±fi]n(±±ffli) 1/2 ' 



The condition for existence of dark solitons is (a + j3i) < 0. In a medium like LiNbO^, a < 0, 
therefore, if \[3\ < \a\ then the dark SP solitons can be observed irrespective of the polarity of external 
bias field. However, photovoltaic constant in some photovoltaic materials, such as, BaTiO^ depends on 
polarization of light 87[ . This means sign of a may be positive or negative depending on polarization of 
light. Therefore, to observe dark SP solitons, polarization of light and external bias must be appropriate 
so that (a + (3) < . It is evident from the expression of dark solitons that, if bias field is much stronger 
in comparison to the photovoltaic field, then these SP dark solitons are just like screening dark solitons. 
On the other hand, if the applied external field is absent, then these dark SP solitons degenerate into 
photovoltaic dark solitons in the closed circuit condition. 

7.8 Self Deflection of Photovoltaic and Screening Photovoltaic Solitons 

In absence of diffusion process, which is usually weak, photovoltaic solitons propagate along a straight 
line keeping their shape unchanged. However, when the spatial width of soliton is small, the diffusion 
effect is significant, which introduces an asymmetric tilt in the light-induced photorefractive waveguide, 
that in turn is expected to affect the propagation characteristics of steady-state photorefractive solitons. 
Several authors 



74 . 17814821 |88[ have examined the self bending phenomenon of PV solitons. The effects of 



81 



82 



higher order space charge field on this self bending phenomenon have been also investigated 
The deflection of PV bright solitons depends on the strength of photovoltaic field. When the PV field E p 
is less than a certain characteristic value, solitons bend opposite to crystal c-axis and absolute value of 



14 



spatial shift due to first order diffusion is always larger in comparison to that due to both first and higher 
order diffusion [79j . When E p is larger than the characteristic value, direction of bending depends on the 
strength of E p and the input intensity of soliton forming optical beam. Self deflection can be completely 
arrested by appropriately selecting E p and intensity of the soliton forming optical beam. 

PV dark solitons experience approximately adiabatic self deflection in the direction of the c-axis of 
the crystal and the spatial shift follows an approximately parabolic trajectory. Nature of self deflection 
of PV dark solitons is different from that of bright solitons in which self deflection occurs in the direc- 
tion opposite to the c-axis of the crystal. Effect of higher order space charge field on self deflection has 
been also investigated [88| . which indicates a considerable increase in the self deflection of dark solitons, 
especially under the high PV field. Thus, the spatial shift due to both first and higher order diffusion is 
larger in comparison to that when first order diffusion is present alone. 

The self deflection of screening PV solitons depends on both the bias and PV fields [8l|. When the 
bias field is positive and the PV field is negative, the screening PV bright solitons always bend in the 
direction opposite to the crystal c-axis, and the absolute value of the spatial shift due to first order 
diffusion term alone is always smaller than that due to both first and higher order. When PV field is 
positive and bias field is negative or both are positive, then the bending direction depends both on the 
strength of two fields and on the intensity of the optical beam. Bending can be completely compensated 
for appropriate polarity of the two fields and optical beam intensity. 



8 Two-Photon Photorefractive Phenomenon 

Three types of steady state photorefractive spatial solitons, as elucidated earlier, owe their existence on 
the single photon photorefractive phenomenon. Recently, a new kind of photorefractive solitons has been 
proposed in which the soliton formation mechanism relies on two-photon photorefractive phenomenon. It 
is understood that the two-photon process can significantly enhance the photorefractive phenomenon. A 



new model has been introduced by Castro-Camus and Magana [91| to investigate two-photon photorefrac- 
tive phenomenon. This model includes a valance band (VB), a conduction band(CB) and an intermediate 
allowed level(IL). A gating beam with photon energy Tiu\ is used to maintain a fixed quantity of excited 
electrons from the valance band to the intermediate level, which are then excited to the conduction band 
by the signal beam with photon energy hu>2- The signal beam induces a charge distribution that is 
identical to its intensity distribution, which in turn gives rise to a nonlinear change of refractive index 
through space charge field. Very recently, based on Castro-Camus and Magana's model, Hou et. al., 
predicted that two-photon screening solitons (TPSS) can be created in a biased nonphotovoltaic photore- 
fractive crystal [q2| and two-photon photovoltaic (TPPV) solitons can be also created in a PV crystal 
under open-circuit condition 93[. Recently, the effect of external electric field on screening photovoltaic 



solitons due to two-photon photorefractive phenomenon has been also investigated 



94, 



8.1 Two-Photon Photorefractive Nonlinearity 

In order to estimate the optical nonlinearity arising out in two photon-photorefractive media, we consider 
an optical configuration whose schematic diagram is shown in figure (8). 

Insert Figure (8) here 
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The electrical circuit consist of a crystal ( could be made of photovoltaic-photorefractive or non 
photovoltaic-photorefractive material), external electric field bias voltage V and external resistance R. 
Vq and Eq respectively denote the potential and electric field strength between the crystal electrodes 
which are separated by a distance d. Assuming the spatial extent of optical wave is much less than d, we 
have Eq — Vq/ d; additionally V/d = E a . 

Therefore, 

V = E a d + JSR = E a d, (46) 



where S is the surface area of the crystal electrodes and J is the current density. The soliton forming 
optical beam with intensity I 2 propagates along the z direction of the crystal and is permitted to diffract 
only along the x direction. The optical beam is polarized along the x axis which is also the direction of 
crystal c-axis and the external bias field is also directed along this direction. The crystal is illuminated 
with a gating beam of constant intensity I\ . The space charge field E sc due to two-photon photorefractive 
phenomenon can be obtained from the set of rate, current, and Poisson's equations proposed by Castro- 
Camus et. al. [91]. In the steady state, these equations are 

fall + Pi){N D - N+) - 7l iV,:iV+ - j R NN+ = 0, (47) 



( Sl h + ^)(N D - N+) + j 2 N(N it - Ni) - 71 - (s 2 I 2 + fa)Ni = 0, (48) 

1 d I 

(s 2 J a + fo)Ni + ~q£- - l2 N(N lt - N t ) = 0, (49) 

e e r -^ = e(N+-N-Ni-N A ), (50) 

J = efiNE sc + k p s 2 (N d - N+)I 2 + eD^-, (51) 
dJ 

—— = 0, J = constant, (52) 
ox 

where Nd , N^, , Na and are the donor density, ionized donor density, acceptor or trap density 
and density of electrons in the conduction band, respectively. Ni is the density of electrons in the 
intermediate state, Na is the density of traps in the intermediate state. k p , fi and e are respectively 
the photovoltaic constant, electron mobility and electronic charge; is the recombination factor of 
the conduction to valence band transition, 71 is the recombination factor for intermediate allowed level 
to valence band transition, 72 is the recombination factor of the conduction band to intermediate level 
transition; f3\ and j3 2 are respectively the thermo-ionization probability constant for transitions from 
valence band to intermediate level and intermediate level to conduction band; s\ and s 2 are photo- 
excitation crosses. D is the diffusion coefficient and I 2 is the intensity of the soliton forming beam. We 
adopt the usual approximations N^ ~ Na and N t — N << Na ■ In addition, we also assume that, 
the power density is uniform at large distance from the center of the soliton forming beam, thus, at 
x — > ±00, I 2 (x — > ±00, z) = constant = I 2oo . Obviously, the space charge field in the above remote 
region is also uniform, i.e., E sc (x — > ±00, z) — Eq. The space charge field can be obtained using standard 



procedure [94|, which turns out to be 
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{hoo+hd){h+hd + ^) 
^ gEa {h+h d ){hoo+h d + ^) 

„ s 2 {ghoo-h){h + hd + 11 ^ 
+ E„- 



{h+hd){sih+Pi) 

_D 7l iV A dh 

MS2 {h + J M )(/aoo + hd + ^) dx ' 



(53) 



where E p — k p Na!r/ e.P' is the photovoltaic field, lid — /W s 2 is the so called dark irradiance. 
g = 1/(1 + q), q = eAljv j SR , Noo = N(x — > ±00). In general the parameter g is a positive quantity and is 
bounded between < g < 1. Under short circuit condition i? = and 9=1, implying that the external 
electric field is totally applied to the crystal. For open circuit condition R — > 00, thus, 5 = i.e., no bias 
field is applied to the crystal. 



8.2 Nonlinear Evolution Equation 

As usual, the optical field of the incident soliton forming beam is taken as ~E = ~x*$(x, z) exp[i(kz — 
cot)], where k = kon e , ko = 2tt/Xo, Xq is the free space wavelength of the optical field, n e is the unperturbed 
extraordinary index of refraction and $ is the slowly varying envelope of the optical field. Employing the 
space charge field as given by equation(53) and following the procedure employed earlier, the nonlinear 
Schrodinger equation for the normalized envelope A(s, £) can be obtained as [94| 

.dA Id 2 A {l + P ){l + a+ \A\ 2 )A 
l dZ + 2ds 2 19 (I + \A\ 2 )(l + a + p) 



ar\- 



(gp- \ A \ 2 )(i 



\A\ 2 )A 



l + \A\ 2 
SaA 



o\a\< 



(1 + \A\ 2 )(1 



a 



\A\ 2 ) ds 



= 0, 



(54) 



where p = hoo/hd, A = y / 2)) "} 2d $, £ = z/(k Q n e xl), s = x/x , P = (k x Q ) 2 (n^r 33 /2)E a , a = 
(k x ) 2 (nir 33 /2)E p , S = (k x ) 2 (njr 33 /2)D/M, r, = fa/{sih + ft), a = = ^ and m = 

\J Mo/eo- ^33 is the electro-optic coefficient of the two-photon photorefractive crystal. Equation(54) can 
be employed to investigate screening, photovoltaic and screening photovoltaic solitons under appropriate 
experimental configuration. 



8.3 Screening Solitons 

For screening solitons, the crystal should be nonphotovoltaic-photorefractive(i.e., a — 0). Assuming that 
the external bias field is totally applied to the crystal(R=0), thus, E a = E and g = 1. Neglecting 



diffusion, the expression for space charge field E sc turns out to be 9ll. |92| 



(hoc+hd)(h+hd+ 11 £ A ) , N 

E sc = E - s ^r-. 55 

{h+hd){hoo+hd + : ^) 

Note that though the gating beam of constant intensity I\ is required to maintain a quantity of excited 
electrons density in the intermediate allowed level, it does not appear in the expression of space charge 
field. The relevant modified nonlinear Schrodinger equation is obtained as 

M 1 ld2A M + P)(l + - + \A\ 2 )A _ 
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Fundamental properties of dark, bright a nd grey solitary wave solutions of this equation have been in- 
vestigated extensively by several authors 93. 196l498l | . Self deflection of these solitons due to diffusion 92 1 
and effects of higher order space charge field on self deflection have been also examined [97[. Jiang et. 
al. have examined temporal behavior of these solitons. They predicted that, in the low amplitude 
regime, FWHM of solitons will decrease monotonically to a minimum steady state value, and that the 
transition time of such solitons should be independent of j3 or soliton intensity and is close to lOT^, where 
Td is the dielectric relaxation time. They also predicted that the temporal properties of dark solitons are 
similar to those of the bright solitons. 



Intensity of bright solitons vanishes at infinity, thus, /200 = p = 0. Therefore, 

.dA Id 2 A (3 



<9£ 2 ds 2 



1 



l + \Af 



A = 0. 



(57) 



Assume a bright soliton of the form A — p 1 ^ 2 y(s) exp(w£), where v is the nonlinear shift in propagation 
constant and y(s) is a normalized real function, which is bounded as < y(s) < 1. The parameter p 
stands for the ratio of the peak intensity of the soliton to dark irradiance l2d- The profile of the soliton 
turns [92J out to be 



s = ± 



r WE 1 

L (l + cr) J 



-l/2pl/2 



,„ [ln(l+pP)-Pln(l+p)]i/ 
while the expression for nonlinear phase shift is given by 



1 



[1 + - ln(l 

P 



P)} 



(58) 



(59) 



From equation(58), it can be easily shown that the bright soliton requires (3 > i.e., Eq > 0. 
Therefore, screening bright spatial solitons can be formed in a two-photon photorefractive medium only 
when external bias field is applied in the same direction with respect to the optical c-axis. FWHM of 
these spatial solitons is inversely proportional to the square root of the absolute value of the external bias 
field. In the low amplitude limit i.e., when \A\ 2 << 1, equation(57) reduces to 



dA Id 2 A 



<9£ 2 ds 2 1 



(1 + a - a\A\ 2 ) A = 0. 



(60) 



Above equation can be exactly solved analytically and the one-soliton solution is 
A(s,£) = P 1/2 sech 



( Ppa\ 1/2 ' 






exp 



p(pa - 2a - 2) ' 
l — ; i 



2(1 + a) 



The spatial width (tfwhm) of these solitons are TpwHM = 21n(l + V2) ( 



1/2 



(61) 



To obtain dark soliton solution, we take A(s, £) = p 1 / 2 y{s) exp[ip£], where y(s) is a normalized odd 
function of s and satisfies the following boundary conditions: y(s = 0) = 0, y(s — > ±00) = 1, and all 
the derivatives of y(s) vanish at infinity. The profile y(s) of these solitons can be obtained [92j using the 
following relationship 
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(62) 



and the nonlinear phase shift /i = — /3. The dark solitons require /3 < 0. In the low amplitude limit 
i.e., when \A\ 2 « 1, cquation(56) reduces to 



.OA Id 2 A 13(1 + p) 



d£ 2ds 2 (1 + a + p) 
The dark soliton solution of above equation turns out to be 



(l + <r -a\A\ 2 )A = Q. 



A{s,Ci = p 1/2 tanh 



ppa 



1 



1/2 



exp 



ip(\+p){pa 



i 



(63) 



(64) 



1/2 



The spatial width ( Tfwhm) of these solitons are tfwhm = 21n(l + V2) (± 
doped strontium barium niobate (SBN) could be a good candidate for observing these solitons, since, 
they have an intermediate level for the two step excitation. In addition to bright and dark solitons, 
equation(56) also predict steady state gray solitons under appropriate bias condition. These screening 
gray solitons were investigated by Zhang et. al. Q|- Properties of these gray solitons are similar to 
fundamental properties of one-photon photorefractive gray spatial solitons. For example, they require 
bias field in opposite to the optical c-axis of the medium and their FWHM is inversely proportional to the 
square root of the absolute value of the bias field. The main difference between one-photon and two-photon 
gray solitons is that one-photon gray solitons rely on one-photon photorefractive phenomenon to set up 
space charge field, while the two-photon photorefractive gray solitons rely on two-photon photorefractive 
phenomenon to set up space charge field. 



8.4 Photovoltaic Solitons 

We consider a photovoltaic photorefractive crystal under open circuit condition (R = 0, g = 0), the 
expression for space charge field from equation(53) reduces to 



E,, 



-E, 



S2h{h + hd 



7i Na ' 
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7i N A 



dh 



[h + hd){hsi + Pi) ps 2 {h + hd) (7 2qo + hd + 2lMA 



dx 



(65) 



Upon neglecting diffusion, the nonlinear Schrodinger equation for the normalized envelope A(s, £) can 
be obtained as 



.OA Id 2 A (l + a+\A\ 2 )\A\ 2 A n 
i— + + arj- ) tt^ A = 0. 



(66) 



d£, 2ds 2 1 (1 + \A\ 2 ] 

Bright, dark and gray photovoltaic solitons of equation(66) have been examined by several au thors [9 
10o| | . Deflection of these solitons and higher order effects have been also given adequate attention 



For bright solitons, we consider a si milar profile as it was considered for screening solitons, and, the profile 

9^. [loot 



of such solitons turns out to be 



s = ±{arj)- 1 ' 2 [ { — [Hi 
J y { P 



+ py 2 ) - f ln(l + p)\ + P y 2 (l - py 2 ) 



-1/2 



(67) 
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where the nonlinear phase shift is given by 



1 OlTjp 

v = arjall ln(l + p)\ H . 

P 2 

The bright soliton requires a > 0, i.e., E p > 0, thus the photovoltaic field should be in the same direc- 
tion with respect to the optical c-axis of the medium. In the low amplitude limit (| A\ 2 << 1), FWHM of 
these solitons are inversely proportional to the square root of the absolute value of the photovoltaic field. 



The profile of dark photovoltaic solitons, following a similar procedure, turns out to be [9 

-1-1/2 



s = ±(—ar/) 



-1/2 



and the nonlinear phase shift is given by 



fi = otr\p 1 + 



1 + P 



(69) 



(70) 



The two-photon photovoltaic solitons require a separate gating beam to produce a quantity of excited 
electrons from the valance band to the intermediate level of the material. Without the gating beam, the 
signal beam cannot evolve into a spatial soliton. By adjusting the gating beam, one can control the width 
as well as formation of two-photon photovoltaic solitons. 



8.5 Screening Photovoltaic Solitons 

Steady state screening photovoltaic solitons are obtainable when an electric field is applied to a photo- 
voltaic photorefractive crystal. These SP solitons result from both the photovoltaic effect and spatially 
nonuniform screening of the applied field. Recently, bright and dark screening photovoltaic (SP) solitons 
have been investigated by Zhang and Liu 94[. The normalized bright field profile y(s) of these solitons 
can be determined ( with 5 — 0) from equation(54), which turns out to be [941 ] : 



s = ±f { Jf° [ln(l + py 2 ) - y 2 ln(l + p)} 

Jy Rl+O") 

+ ^[l n (l +pf) - y 2 ln(I +p)} + aripy 2 (l - y 2 )y 1/2 dy. (71) 

The nonlinear phase shift v of these solitons is given by 

v = + - ln(I + p)] + a W [l - - ln(l + p)} + (72) 

1 + a p p 2 

Please note that unlike bright screening solitons, in the present case, it is not necessary that the value 
of (3 should be positive. However, the sign of a and (3 should be such that the curly bracketed term in 
equation(71) is positive. From equation(54) the normalized dark field profile y(s) can be obtained as 94|: 
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Jy (1 + P + 0") 

1 + p f l + py 2 \ r 2g(l +gp) 2 

2.(i + 5 p) ln /i^\^ 2 _ 1)2]rl/2 ^ (?3) 



p V i + P 

where the nonlinear phase shift is given by 

(i = -g/3 + a VP (l - g) (l + . (74) 

By setting a — and g = 1 in cquation(71), we recover bright solitons of equation(58). Similarly, by 
setting a — and g = 1 in cquation(73), we recover dark solitons of equation(62). In addition, by 
taking f3 = and g = in equation(71) i.e., in open circuit realization, we recover bright PV solitons of 
equation(67). Similarly, by setting g = and j3 — 0, we recover dark PV solitons of equation(69) from 
equation(73). 

As pointed out by Zhang and Liu [94], these two-photon SP solitons may be considered as the unity 
form of two-photon screening and two-photon photovoltaic solitons under open circuit realization. If the 
biased field is much stronger in comparison to the photovoltaic field, then the screening photovoltaic 
solitons are just like screening solitons. If the applied field is absent, the screening photovoltaic solitons 
degenerate into the photovoltaic solitons in the open circuit condition(/3 = 0,g — 0). In other words, 
the open circuit photovoltaic solitons or screening solitons are special cases of the screening photovoltaic 
solitons. Equation(73) also predicts the existence of two-photon photovoltaic solitons when 5 = 1 and 
f3 = i.e., two-photon photovoltaic solitons in closed circuit realization. 

Before closing this section, a brief comment on gray two-photon screening PV solitons in biased two 
photon phototvoltaic crystals seems inevitable. Equation(54) predicts the existence of such solitons 
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The properties of these gray solitons, such as, their normalized intensity profiles, intensity F WHM , 
transverse velocity and transverse phase profiles have been discussed in detail by Zhang et. al. 1021 ] . 



They become narrower as the grayness parameter m decreases for a given normalized intensity ratio p. 
However, the soliton width generally decreases and transverse velocity generally increases with intensity 
ratio p. In addition, soliton phase varies in a very involved fashion across transverse direction and the 
total phase jump of these solitons exceeds 7r for relatively low value of the grayness parameter m. 



9 Vector Solitons 

Thus far, we have discussed optical spatial solitons which are solutions of a single NLS equation. These 
solutions are due to a single optical beam with a specific polarization and the polarization is maintained 
during propagation. However, always this specific picture may not hold good. Two or more optical 
beams may be mutually trapped and depend on each other in such a way that each of them propagates 
undistorted. Thus, several field components at different or same frequencies or polarizations may interact 
and yield shape preserving propagation. In order to discuss such cases, we need to solve a set of coupled 
NLS equations. Shape preserving solutions of this set of coupled NLS equations are called vector solitons. 
Only in specific cases, the constituents of these solitons are vector fields associated with solitons. In 
general, they are multi component in nature. 
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9.1 Two Component Incoherently Coupled Vector Solitons 

Among spatial solitons interaction, pairing of two spatial solitons has been always an intriguing and ex- 
tensively investigated issue. When two such soliton forming beams propagate, they interact through cross 
phase modulation (XPM) and induce a refractive-index modulation created by both beams. Two beams 
are mutually trapped and depend on eac h other in such a way that each of them propagates undistorted. 



Very recently, vector screening solitons 104| - |l07j | have been investigated, that involve two polarization 
components of an optical beam which are orthogonal to each other. Depending on the symmetry class of 
the crystal and its orientation, these solitary beams obey cross or self coupled vector systems of dynamical 
equations. 

A new t ype of steady state incoherently coupled soliton pair was discovered in biased photorefractive 



crystals 1081 ] . which exists only whe n the tw o soliton forming beams possess same polarization and fre- 



quency and are mutually incoherent 108 123] . These solitons can propagate in bright-bright, dark-dark, 
bright-dark and gray-gray configurations, and they can be realized in simple experimental arrangement 
with two mutually incoherent collinearly propagating optical beams. Since two beams are mutually inco- 
herent, no phase matching is required and they experience equal effective electro-optic coefficients. The 
idea of two incoherently coupled solitons has been generalized and extended to soliton families where 
number of constituent solitons are more than two. Such incoherently coupled families can be established 
provid ed they have sa me polari zation, wavelength a nd a re mutually incoherent. Bright-bright, dark- 
dark 118j . bright-dark 119Ml2l| as well as gray-gray 113j configurations have been investigated. These 
multi component solitons are stable. In next few sections, we confine our discussion on two component 
spatial photorefractive vector solitons which are co-propagating and overlapping. 

9.1.1 Coupled Solitary Wave Equations Due to Single Photon Phenomenon 

To start with, we consider a pair of optical beams which are propagating in a photorefractive crystal ( 
the crystal could be PV-PR or non PV-PR) along z-direction. They are of same frequency and mutually 
incoherent. The optical c-axis of the crystal is oriented along the x direction. The polarization of both 
beams is assumed to be parallel to the x-axis. These two optical beams are allowed to diffract only along 
the x-direction and y-dynamics has been implicitly omitted in the analysis. For the sake of simplicity 
the photorefractive material is assumed to be lossless. The perturbed refractive index along the x-axis is 
given by — n\ — n^r^E sc . The optical fields are expressed in the form E\ = x $1(2;, z) exp(ifcz) and 
E2 = 1$<&2(x, z) exp(ifcz), where $1 and $2 are slowly varying envelopes of two optical fields, respectively. 
It can be readily shown that the slowly varying envelopes of two interacting spatial solitons inside the 
photovoltaic PR crystal are governed by the following evolution equations: 

For relatively broad optical beams and under strong bias condition, the space charge field can be 
obtained from cquation(ll) as 

E sc = E I -f^+E P I -^I, (77) 

where in the present case / = I(x, z) is total power density of two optical beams, 1^ is the total power 
density of soliton pair at a distance far away from the center of the crystal i.e., 1^ = I(x — > ±00). E 
is the value of the space charge field at far away from the beam center i.e., x — > ±00. For two mutually 
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incoherent beams, total optical power density I can be written as I(x,z) = n e /(2fyo)(| < I > i | 2 + l^l 2 ) i-e., 
sum of Poynting fluxes. Substituting the expression of E sc in equations(75) and (76), we derive the 
following dimensionless dynamical equations for two soliton forming optical beams: 



OA, 1 d 2 A, 

— - H 

«9£ 2 ds 2 
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(78) 



where /3 = (k x ) 2 (n^r 33 /2)E , a = (k Q x ) 2 (n^r 33 /2)E p , Aj = ^ 



- y 2^®i' £> s and /° are defined 
earlier. Above set of two coupled Schrodinger equations can be examined for bright-bright, brig ht-dark, 
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dark-dark, gray-gray screening, photovoltaic as well as screening photovoltaic solitons 
be easily shown that when the total intensity of the two coupled solitons is much lower than the effective 
dark irradiance, the coupled soliton equations reduce to Manakov equations. The dark-dark, bright-bright 
and dark-bright soliton pair solutions of these Manakov equations can be obtained under appropriate bias 
and photovoltaic fields [117| . 

With the growing applications of self focusing and spatial s olitons in modern technology, several mathe- 
matical methods have evolved to address soliton dynamics 124 - 129^ . In particular, Christodoulides et. 
al. 108] have developed a very efficient method to numerically solve a set of coupled equations which has 



been employed extensively to investigate coupled solitons in PR media. In this method, two coupled equa- 
tions are converted to one ordinary differential equation which is then numerically solved to obtain soliton 
profiles. The main difficulty with this method is its inability to capture the existence of a large family 
of stable stationary solitons. We will discuss more on this in the latter part of the article. What follows 
in the next section is a discussion on incoherently coupled solitons employing the method of reference 108. 



9.2 Incoherently Coupled Screening Vector Solitons 

In this section we consider incoherently coupled bright-bright, dark-dark as well as bright-dark screening 
solitons in a biased nonphotovoltaic photorefractive crystal. These coupled solitons were identified by 
Christodoulides et. al. [108]- The parameter a = 0, since, the crystal is nonphotovoltaic, hence, relevant 
coupled Schrodinger equations are as follows 



'% + ^-« 1+ ' WiVW) =0 ' (79) 

We first consider a bright-bright soliton pair for which 1^ = p = 0. Expressing stationary soliton 
solutions of the form A\ = p 1 ^ 2 y(s) cos 9 exp(«/i^) and A2 — p 1 ^ 2 y(s) sin# exp(i/i£), where fj, represents 
nonlinear shift of the propagation constant, and y(s) is a normalized real function between < y(s) < 1. 
The parameter 9 is an arbitrary projection angle which ultimately decides the relative power of two 
components. Substituting A\ and Ai in equation(79), we get the following ordinary differential equation 



i? 

ds 



2/?: 



0. 



(80) 



and 



M =-(-)ln(l+p). 
P 



(81) 
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Earlier in section (7.3), it was shown that, above equation admits bright solitons when /3 i.e., Eq is 
positive. The, same condition holds good for bright-bright pair which can be obtained by numerically 
solving equation (80). A typical bright-bright pair has been depicted in figure (9). 

Insert Figure (9) here 

For dark-dark pairs, 1^ and p are finite. We express A\ and A 2 as A\ = p 1 ^ 2 y(s) cos9 exp(ip^) and 
A2 = p 1 ^ 2 y(s) sin9 exp(ip{;), with \y(s)\ < 1. Thus, we have 

,2 

-1 - 2f xy -2/3(1 + p)— V — = 0, (82) 
ds z 1 + py z 

and 

p = (83) 

Equation (82) can be solved for dark-dark soliton pairs provided /3 i.e., Eq is negative. A typical 
dark-dark pair has been depicted in figure (10). 

Insert Figure (10) here 

For bright-dark soliton pair, we express A\ and A2 as A\ — p x l 2 f(s) exp(i/z£) and A2 = p 1 ^ 2 g{s) exp(w£), 
where f(s) and g(s) respectively represents envelope of bright and dark beams. Two positive quantities 
p and p represent the ratios of their maximum power density with respect to the dark irradiance Id- 
Therefore, bright-dark soliton pair obeys following coupled ordinary differential equations: 



and 



ds 2 



fg 

ds 2 



- 2 



0(1 + P) 



1+pf 2 



pg 



- 2 



P0- + P) 



1 + pf 2 



pg 



/ = o, 



g = 0. 



A particular solution of above equations can be obtained using the simplification f 2 



(84) 



(85) 



1. Em- 



ploying appropriate boundary conditions, the nonlinear phase shifts p and v are obtained as: p = 
— ^ ln(l + A) and v = —/3, where A = (p — p)/(l + p). Whe n peak int ensities of two solitons are approx- 



imately equal, A << 1 
Ax = P 1 / 2 sech[((3A)^ 2 s} exp[ 



The approximate soliton solution [1051 Il08| in this particular case is given by 
-i/3(l-A/2)£], and A 2 = p 1/2 tanli^A) 1 /^] exp[-i/3f]. These two soliton 
solutions are possible only when the product (/3A) is a positive quantity. A typical bright dark soliton 
pair is depicted in figure (11). 



Insert Figure (11) here 



9.2.1 Missing Bright Screening PV Solitons 

In the preceeding section, we have followed the procedure developed in reference 108 and assumed a 
particular type of ansatz for the bright-bright pair. According to this method, the normalized power 
P\ and P2 of two solitons of the bright-bright pair are Pi = psin 2 9 and P2 = pcos 2 #. Therefore, for 
a given Pi, 9 has a fixed value, hence, the power of the other component has only one possible value 
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which is unique. Or in other words a composite soliton can exist only with a single power ratio. In 
this section, we will show that for a given power of on e component, the other component can exist with 
different power. Thus, the method of reference [xOs| fails to identify a large number of bright-bright 
solitons in a two component configuration. Therefore, our goal, in this section, is to demonstrate the 
existence of a new very large family of two-component composite screening photovoltaic spatial solitons 
in biased photovoltaic-photorefractive crystals which were not identified by the method of reference 108. 
For bright solitons, p = 0, hence, relevant coupled equations for bright-bright screening PV solitons in 
biased PV-PR crystals [64j are 

dAj l &Aj [a{\A^ + \A^)~p]A ] 

di + 2 ds 2 + (I + |^f + |^3-, | 2 ) ' ( } 

In order to analyze the behavior of these coupled solitons, we assume solutions of the form 

^M = ^MexpH%£), (87) 
By virtue of use of equation(87) in equation (86), we obtain following equations: 

da, i {on-V i d 2 *, 



dVj dVL 3 dVj 1 d 2 fl 



<9£ ds ds 2 ds 2 



= 0, (89) 



where II = "i+^^J-f j^p^ ■ The last three terms of (88) determine the behavior of the eikonal flj 
i.e., the convergence or divergence of two optical beams. The fourth term in this equation represents 
nonlinear refraction while the third term determines diffraction. Equation (89) determines the evolution 
of the beam envelope ^j. In equation (88), II represents the contribution from nonuniform screening of 
the applied electric field and photovoltaic properties of the crystal as well. 

Lowest order localized bright solitons, for which light is confined in the central region of the soliton, 
obey ^j(s = 0) = $j mol and &j — as s — > ±oo. The fundamental solutions of coupled Schrodinger 
equations, in a self focusing Kerr medium, are represented by sech functions. Equation (86) is a modi- 
fied nonlinear Schrodinger equation(MNLSE) in saturating media. Due to the saturating nature of the 
medium, it is expected that the fundamental soliton solutions will not be exactly sech function. However, 
in many nonlinear optical problems involving NLSE and MNLSE in Kerr, cubic qu intic, nonlocal and 



saturat ing m edia, approximate solutions have been obtained using Gaussian [13CH135} or super Gaussian 
ansatz |135j . The motivation of employing such ansatz is two fold. Firstly, particularly true for Gaussian 
ansatz, mathematical formulations become easy. Secondly, numerically computed exact solutions are not 
widely different from Gaussian profiles in many cases. Thus, though approximate, Gaussian profile still 
provides good approximation to the problem. Hence, the solutions of above equations are taken to be 
Gaussian with amplitude and phase of the following form 



#,-(s,£) = exp 



2r 2 (0/, 2 (0 



(90) 



s 2 



%M = ^(0 + ^(0, ( 91 ) 
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and 



where, ^ j represents the peak power of the component of bright-bright solitons, Tj is a positive 
constant, fj(£) is variable spatial width parameter; rjfj(£) is the spatial width of these solitons and 
4>j(£) is an arbitrary phase function. For a pair of nondiverging solitons at £ = 0, we shou ld ha ve fj = 1 



and -Jjj^ = 0. Substituting for VE^ and Slj in (88), using paraxial ray approximation |86l. Il26l Il27| and 
equating coefficients of s 2 from both sides of (88), we obtain 



Hi _L n[*i ^ ^ 1 



C \^72+ 2 A ~, 72- ( 93 ) 



where Pj = ^ 2 and C — 2 (a + (3). The solutions of above equation should give stationary and 
non-stationary coupled solitons of (86) for given set of power and spatial width. 

9.2.2 Stationary Composite Solitons 

In order to identify stationary composite solitons, we need to locate equilibrium points. The equilibrium 
points of equation(93) can be obtained from following equations 

+ 7. , „V „,2 =0, (94) 



and 



(1 + Pi + P 2 



"7 - c ( 4 + 4 1 2 = °- ( 95 ) 



From above equations, it is obvious that r± = r 2 = r is the condition for existence of stationary 
coupled solitons, where r is a constant. Therefore, composite solitons with different spatial widths 
cannot propagate as a stationary entity. The existence equation of coupled solitons turns out to be 

Or' . <' + * + . (96) 

r\ + f 2 

Obviously, for a bright-bright pair C should be positive. From the existence equation, it is evident 
that though spatial widths of each of the two components are equal, their respective peak power can have 
different value. Equation (96) is equivalent to a quadratic equation in Pi, the root of which is obtained 
as 



Pi = 



{Cr 2 - 2P 2 - 2) ± J(Cr 2 - 2P 2 - 2) 2 - 4(Pf + 2P 2 + 1 - Cr 2 P 2 ) 



12. (97) 



Pi and P 2 are real and positive, hence, always Cr 2 > 4 . For a given value of C, this relationship 
dictates a minimum width for the propagating soliton pair. The variation of Pi with P 2 for different 
values of r has been depicted in figure 12(a)-(b). Each point on any curve of these figures represents a 
stationary composite soliton with a definite spatial width and peak power. 



Insert Figure (12) here 
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An important issue is, whether for a given peak power of one of the component, the other component 
exists with only one or multiple values of peak power. This issue can be settled from figure (13) which 
shows the variation of width with peak power of one component keeping peak power of other component 
constant. It is evident from the figure that Pi has a range of values for a fixed P 2 , thus, with fixed P 2 
other component can exist with different values of P\ . At this stage it is worth pointing out that these 
solitons cannot be identified with the method employed in ref. 108. 

Insert Figure (13) here 



9.2.3 Degenerate Bright Screening PV Bistable Solitons 

We take up a degenerate case in which peak power of two components is same and having same spatial 
width. Setting Pi = P2 = P in (97), we obtain a quadratic equation of P, the solution of which is 
obtained as P = [Cr 2 — 2 ± yj (Cr 2 — 4) 2 — 4]/4 , which implies that spatial width r of each component 
of the composite soliton should be greater than i.e., a two-component composite soliton whose 

individual spatial width is less than above value cannot propagate as a self trapped mode. In figure (14) 
we have displayed variation of r with peak power P. From figure, existence of a bistable regime 



evident i.e., two sets of soliton pairs exist with same spatial width but having different peak power and 
consequently different peak amplitude. Only this degenerate case possesses bistable property. 

Insert Figure (14) here 



9.2.4 Numerical Simulation 

To verify the predictions of foregoing analysis, it is essential to perform numerical simulation. Equation 
(86) has been solved numerically using the split step Fourier beam propagation method 3]. To begin 
with, we look for behavior of the soliton at low power. From figure 12(b), we choose Cr 2 — 4 and r = 1/5 
and select different points from the curve leveled with this value. Chosen points have following values of 
Pi and P 2 , particularly, (i)P x = 0.5, P 2 = 0.5, (ii)Pi = 0.6666, P 2 = 0.3333, (iii)Pi = 0.8333, P 2 = 0.1666 
and (iv)Pi = 0.909, P 2 = 0.0909. It must be emphasized that each point corresponds to a stationary com- 
posite soliton, a paraxial theory prediction. With these parameter values, we launch two Gaussian optical 
beams A\ = \/Pi~exp (^—^pi^j and A 2 — -\/P^exp (—^5^ in (86). Both Gaussian beams acquire solitonic 
shape asymptotically without major modification within very small distance and then they propagate 
almost as a stationary composite soliton. The behavior of two Gaussian spatial solitons corresponding to 
each of these points has been depicted in figure (15). 



Insert Figure (15) here 

It is evident from these figures that a soliton with large power can trap another soliton whose power 
is much lower and both can propagate as a stationary bound state. At this stage it would be appropriate 
to cite one practical example. Consider a BaTiO^ crystal at a wavelength Ao = 0.5/zm with following 
crystal parameters |6 ll : n e = 2.365, r^ = 80 x 10~ 12 m/V, E p = 10 5 V/m. We take arbitrary spatial 
scale Xo = 40/zm and E a — 5.8 x 10 4 V/m. With these values, we find, f3 = 18.4, a = 31.58 and 
C = 2(a + [3) = 99.96. For Cr 2 =4,r« 0.20, thus, in natural unit the intensity FWHM ( i.e.,1.665r) of 
two components is found to be 13.32/irn . An important point to note is that the present investigation 
is also valid for LiNbO^. The LiNbOs parameters could be taken as n e = 2.2, 7-33 = 30 x 10 _12 m/V 
and \E P \ = 4 x 10 6 V/m. However, it should be pointed out that, while E p could be either positive or 
negative for BaTiO^, depending on the polarization of light, the experimental results show that E p is 
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always negative for LiNbO^. Therefore, with proper choice of the value of E a , C = 2(q + j3) can be 
made positive, and hence bright-bright coupled soliton pairs is also observable in LiNbOz- 

We now proceed to obtain bistable composite solitons numerically. From figure (14) we choose two 
points A and B, each corresponds to one pair of composite soliton. These two points have equal values 
of Cr 2 (= 4.2) but two different values of soliton peak power P. We have taken C = 50, hence, r = 0.290 
for both A and B. Peak power P for points A and B are 0.321 and 0.779, respectively. Numerically 
obtained dynamic evolution of a pair of composite bistable solitons corresponding to above values has 
been demonstrated in figure (16). 

Insert Figure (16) here 



In the figure, the upper panel represents one composite soliton while the lower panel represents another. 
From figure, it is evident that both components of the composite in the upper panel remain absolutely 
stationary as they propagate, an example where paraxial theory prediction complies with high accuracy. 
However, in the lower panel, both components of the composite propagate as stable self trapped mode 
though they keep on gentle breathing. Hence, in this case, though prediction of paraxial approximation 
is not very accurate, yet it is able to capture the overall features broadly. Finally, we have found that 
bright-bright pairs are stable against small perturbation in peak amplitude and spatial width. Both 
paraxial theory and numerical simulation show that identified composite solitons are stable. 



9.3 Incoherently Coupled Solitons Due to Two-Photon Photorefractive Phe- 
nomenon 

To investigate coupled solitons in two-photon photorefractive media, the required optical configuration 
is very similar to the one discussed in sec 8.1. The only difference between this case and the earlier 
one is that, in the present case there are two incoherent soliton forming beams whose polarization and 
frequencies are same, whereas in the former case there is only one soliton forming beam. As usual, the 
optical fields are expressed in the form E\ — x^±(x, z) exp(ikz) and E2 = X<f>2(%, z)exp(ikz), where $1 
and $2 are slowly varying envelopes of two optical fields, respectively. The coupled Schrodinger equations 
for the normalized slowly varying envelopes of two optical fields can be described as 



.dU ld 2 U (1 + P ){1 + <j + \U\ 2 + \V\ 2 )U n{gp-\U\ 2 -\V\ 2 ){l + <j+\U\ 2 + \V\ 2 )U 



di 2 ds 2 



(l + a + p){l + \U\ 2 + \V\ 2 ) 



(l + \U\ 2 + \V\ 



= 0, 

(98) 



and 



+ M + p)(i + ° + \u\ 2 + \v\ 2 )v v(gp - P\ 2 ~ M 2 )(i + z + \ u \ 2 + \ v \ 2 ) v = n 

% di^2ds 2 91 (l + a + p)(l + \U\ 2 + \V\ 2 ) " (1 + \U\ 2 + \V\ 2 ) 

(99) 

where U = \J 2 ^i 2d ( ^ >1 ' ^ = J 2-noi 2d ® 2 ' P arame ters p, £, s, f3, a, r), a have been defined earlier. 
Incoherently coupled solitons in two-photon photorefractive media have recently received tr emendou s 



attention, since, the dynamics of these solitons can be controlled by a separate gating beam 137H152j . 
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Properties of these solitons can be investigated using equations (98) and (99). In next few sections we will 
present a brief description of these solitons. 



9.3.1 Incoherently Coupled Two-Photon Photovoltaic solitons Under Open Circuit Con- 
dition 

In this section, we discuss the existence and nonlinear dynamics of two-component incoherently coupled 
composite solitons in two-photon photorefractive materials under open circuit condition. In the steady 
state regime, these incoherently coupled solitons can propagate in bright-dark, bright-bright and dark- 
dark configurations. These photovoltaic soliton families can be established provided that the carrier 
beams share same polarization and wavelength, and numerical simulations show that these solitons are 
stable for small perturbation on amplitude. For photovoltaic solitons under open circuit configuration, 
g = and f) = 0, hence, relevant Schrodinger equations are 

1 + 2d s z +va (i + |c/| 2 + |y| 2 ) ~ u ' 1 UUj 

dv id 2 v (\u\ 2 + \v\ 2 )(i + « + \u\ 2 + W\ 2 )v , lon 



9.3.1.1 Bright-dark solitons 

We first discuss the properties of photovoltaic bright-dark soliton pairs. To obtain the solution for a 
bright-dark soliton pair, the normalized envelopes U and V are expressed as 



U=p 1 ' 2 f{8)exp(irt), 
V = p 1 ^ 2 g{s)exp(iv£) . 



(102) 
(103) 



In above expressions, f(s) and g(s) are real functions, which correspond to the bright and dark pro- 
file, respectively. These real functions are bounded i.e., < f(s) < 1 and < g(s) < 1. Also, p and 
p respectively represents the ratio of solitons maximum intensity to the dark irradiance l2d- Inserting 
expressions (102) and (103) in equations (100) and (101), we obtain 
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(104) 



(105) 



We look for a particular solution which satisfies the condition f 2 + g 2 = 1 . Nonlinear propagation 
constants fj, and v can be determined using appropriate boundary conditions. The value of these turn 
out to be 
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(106) 
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-(B) 

is desirable. In order to do that we integrate equation(104) to obtain 



where A = ( J . At this stage a comment on the sign ol a for the existence of bright-dark solitons 
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d/. (108) 



The sign of the integrand within the third bracket depends on the parameters p, p and a. For a given 
set of experimentally relevant values of the aforesaid parameters, for example, when p = 3.95, p = 4.0 and 
g = 10 6 , the integrand within the third bracket is positive. Thus, for above set of parameters the existence 
of dark-bright solitons requires a < i.e., E p < 0. For illustration, we consider a LiNbO^ crystal with 
the following parameters n e — 2.2 and r^d. = 30 x 10~ 12 mV r_1 at wavelength Ao- Other parameters are 
taken as E p = -4 x 10 6 Vmr 1 , s x = 3 x 10- 4 m 2 W -1 s -1 , 71 = 3.3 x 10- 17 to 3 s~\ N a = 10 22 to~ 3 , j3 x = 
0.5s -1 , fa — 0.5s -1 and s 2 = 3 x 10~ 4 m, 2 W~ 1 s~ x . The gating beam intensity I± = 10 6 W/m 2 , the 
scaling parameter xq — 0.5/xm, therefore, a w —22.2 and 77 = 1.67 x 10~ 4 . The value of tr can be 



controlled by modulating the dark irradiance artificially using incoherent illumination 93j and for the 
present investigation we take a = 10 6 and p — 4. A typical bright-dark pair has been depicted in figure 
(17). 

Insert Figure (17) here 

In order to examine the influence of the gating beam on these solitons, we have numerically computed 
profiles of these solitons at two different values of the gating beam intensities. This has been depicted in 
figure (18). With the change in I\ , the width of each component changes. 

Insert Figure (18) here 

9.3.1.2 Bright-bright solitons 

We now investigate two component bright-bright solitons. In this case, intensities of both soliton forming 
optical beams vanish at infinity i.e., as s — > ±00, /200 = 0. The soliton solution is now expressed in terms 
of normalized envelopes U and V as 

U =p 1/2 y(s)cos9exp{int), (109) 

V = P 1/2 y(s)sin8 exp(i/<), (110) 

where p represents the ratio of the peak intensity to the dark irradiance / 2 rf, P- is the nonlinear shift of 
the propagation constant, y(s) is the normalized real function which is bounded as < y(s) < 1 , 6 
is an arbitrary projection angle which describes relative strength of two components of the composite. 
Substitution of expressions (109) and (110) in either of equations (100) or (101) yields the following 
differential equation, 



il=2(p- a W )y - 2a V py 3 + = 0. ! I I I i 

as 1 1 + py A 

Integrating above equation once, we obtain, 



f s J = 2(p - a W )(y 2 - l)-anp(y*-l) + ^ log l^-^- ) . CI.] 2) 
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Making use of the boundary conditions y(±oo) = and y(±oo) = 0, we can easily obtain p as 



p = ar\a 

Inserting equation (113) in (112) we get, 



' _ log(l + p) 



^ (H3) 



= a{^[log(l +py 2 ) — y 2 log(l + p)} + VP y 2 (l -y 2 )}. (114) 
as J p 

From equation (114) we can easily show that the quantity within the curly bracket in the right hand 
side is positive for all the values of y 2 (s) i.e., < y(s) < 1, therefore, we easily conclude that a > for 
bright-bright solitons. In order to investigate bright-bright soliton pair, we take a Cu:KNSBN crystal, 
whose parameters at A = 0.5pm are taken as n e = 2.2, r 33 = 200 x lO^^mV^ 1 and E p = 2.8 x IQ^Vm^ 1 . 
The scaling parameter x = 10pm, p — 10 and 9 = 30°. Other parameters for the bright-bright soliton 
configuration are: a — 22.2, r\ ~ 1.5 x 10~ 4 and a = 10 4 . Figure (19) depicts the normalized intensity 
profile of the photovoltaic bright-bright soliton pair. 

Insert Figure (19) here 

9.3.1.3 Dark-dark solitons 

Properties of dark-dark soliton pairs can be analyzed following similar procedure as elucidated in previous 
sections. In the case of dark type profiles, there is a constant intensity background i.e., 7 2 oo ^ 0, therefore 
we express normalized envelopes U and Va,s 

U = p 1/2 y(s)cos6e?ip(i(iZ), (115) 
V = p 1/2 y{s)sin6cxv{i^), (116) 



where, p = hoo/hd and < y(s) < 1. As usual p is the nonlinear shift of the propagation constant, 
and 9 is the projection angle. Furthermore, inserting equations (115) and (116) in either of equation 
(100) or (101), we obtain 

g = 2(p a W )y 2ar 1 py 3 + (117) 

Above equation can be solved easily adopting numerical procedure after evaluationg the nonlinear 
propagation constant p. An important point to note is that dark-dark soliton pairs require a < 0. For 
illustration, we take a LiNbOz crystal with n e = 2.2 and r 33 = 30 x lO^^mV^ 1 at wavelength A = 
0.5^m. Other parameters are E p = -4 x VPVmr 1 , s 1= 3x lO-^W' 1 s' 1 , 71 = 3.3 x 10" 17 m 3 S - 1 , 
N A = 10 22 TO~ 3 , /3i = 0.5s _1 ,/3 2 = 0.5s" 1 and s 2 = 3 x 10 _4 m 2 VF _1 s _1 . The gating beam intensity Iiis 
taken to be l0 6 W/m 2 and the scaling parameter xo = 0.5pm. Therefore, a w —22.2, 77 = 1.67 x 10~ 4 . 
We take a = 10 6 and p — 4. A dark-dark soliton pair is depicted in figure (20). Numerical simulation 
confirms that these solitons are robust, do not break up or disintegrate if small perturbation in amplitude 
is introduced. 

Insert Figure (20) here 



31 



10 Conclusion 



We have presented a brief review of the recent developments in the field of optical spatial solitons in 
photorefractive media. In relatively short time this topic has achieved tremendous success in theory as 
well as in experiments. We have considered fundamental properties of three types of solitons, particularly, 
screening, photovoltaic and screening photovoltaic solitons and described different methods to investigate 
them. For each type of soliton, three different configurations i.e., bright, dark and gray varieties have 
been considered. Self bending of these solitons due to diffusion and effect of higher order diffusion on 
self bending phenomenon are also highlighted. Besides single photon photorefractive phenomenon, lately 
the two-photon photorefractive phenomenon has become a topic of intense research since the PR effect 
can be controlled with a separate gating beam. Mechanisms of formation of PR solitons due to single 
photon as well as two-photon photorefractive processes have been discussed. Interaction of solitons is an 
extremely important topic which could be exploited to fabricate all optical switching devices. We have 
discussed important properties associated with interaction of these solitons. Vector solitons, particularly, 
incoherently coupled solitons due to single photon and two-photon photorefractive phenomena have been 
highlighted. Existence of some missing solitons pointed out. Several properties discovered so far for these 
solitons are universal and applicable to other branches of solitons. 
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Figure 1: Band transport model 




Figure 2: Normalized intensity profile of bright spatial solitons for j3 = 43, xq — 20 //m;p = 0.1 and 10. 
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Figure 3: Variation of spatial width (tfwhm) of solitons with power p. Figure shows existence of bistable 
solitons. 




Figure 4: Normalized intensity profile of a pair of bistable solitons. Both solitons have same spatial width 
but they possess different peak power. 
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Figure 5: Normalised intensity profile of dark spatial solitons for j3 = —43, xq = 20^m; p = 1 and 2. 




Figure 6: Existence curve of stationary solitons for different /3. 
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Figure 7: Spatial shift of soliton centre as it propagates through the photovoltaic crystal. Aqq = 0.42 
and r = 0.298. 
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Figure 13: Variation of peak power Pi of one of the component of the composite soliton with spatial 
width r while the peak power of other component P2 is constant. 
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Figure 14: Variation of peak power P of the degenerate composite soliton with spatial width r. Nature 
of the curve signifies existence of bistable property of solitons. 
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Figure 15: Stable propagation of composite solitons as obtained by direct numerical simulation. Cr 2 = 4, 
r = 1/5. (a) Pi = 0.5, P 2 = 0.5, (b) P 1 = 0.6666, P 2 = 0.3333, (c) Pi = 0.8333, P 2 = 0.1666, and (d) 
Pi = 0.909, P 2 = 0.909 . Left panel |Ai| 2 and right panel \A 2 \ 2 . 
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Figure 16: Propagation of bistable composite solitons. Peak power and width of these solitons have been 
chosen from point A and B of figure (14) which corresponds to Cr 2 = 4. Upper panel corresponds to 
point A and lower panel corresponds to point B. Solitons of both upper and lower panels have same width 
i.e., each component has spatial width r — 0.290. Peak power of each component in the upper panel 
P = 0.321. Peak power of each component in the lower panel P = 0.779. 
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Figure 17: Soliton components \U\ 2 and \V\ 2 of the bright-dark soliton pair for p = 4,p — 3.95, a = 10 6 
and a = —22.2. Gating beam intensity I\ = 10 6 W/m 2 and calculated value of A = —0.01. 




Figure 18: Soliton components \U\ 2 and \V\ 2 of the bright-dark soliton pair for two different values of 
the gating beam intensity I\. Values of different parameters are p = A,p = 3.95, a = 10 6 and a = —22.2. 
Solid line for h = 10 5 VF/m 2 , dashed line for h = 10 6 W/m 2 . 
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Figure 19: Soliton components \U\ 2 and \V\ 2 of the bright-bright soliton pair when p — 10 and 9 = 30°. 
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